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ABSTRACT 


The factors affecting the stability and the maneuvera- 
bility of a dynamically unstable Ship were studied using 
the linear and non-linear equations of the motion of the 
ship. The DSL/360 Language was used to simulate the dynam- 
ics and to study both unstable motions and standard 


Mmaxarecuvers. 
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Note: 


1) Symbol of prime ('), except where otherwise stated 
indicates the values of an item non-dimensionalized. 


2) symbol of dot (°) indicates differential value of 
an item with respect to time. 


oo oubscript f indicates the value of an item of fin. 


4) Subscript h indicates the value of an item of ship's 
nude. 


a subscript ry indicates the value of an item of the 


large deadwood included rudder. 


6) Subscript fr, indicates the value of an item of the 


small deadwood or rudder only. 
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te UNGAR ODUCTION 


A ship is said to be dynamically stable on a straight 
@eourse Or in a turn of constant curvature if, when slightly 
disturbed from it's steady motion, it will soon resume that 
same motion along a slightly shifted path without any cor- 
recting control being applied. A ship which is dynamically 
unstable in straight line motion cannot maintain straight 
line motion when there is no control. However, with regard 
to maneuverability, a ship which is too stable will not 
turn as tightly as a somewhat less stable ship so that a 
highly stable ship may compromise maneuverability. There- 
fore, a ship should be designed for a moderate amount of 
stability so as to be able to go straight without compro- 
mising maneuverability too much. 

With these factors in mind then, the stability of five 
different ships was studied in Part III and maneuverability 
Miers lV¥. The first objective of thas thesis was to 
investigate the stability of a ship with movable fin. The 
second objective was to investigate the maneuverability of 
an unStable ship in standard maneuvers to compare with those 
Sia Stable ship. A third objective was to apply DSL/360 


Language to simulate the ship motions. 
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II. EQUATION OF SHIP MOTION IN THE HORIZONTAL PLANE 


As developed in Ref. 2, the general motion of a ship is 


feat Of a rigid body in six degrees of 
sravity and bucyancy forces as well as 
and hydrodynamic reactions and subject 
Sener disturbances or excitations. 
Considering axes fixed to the body 


the principal axes of inertia but with 


freedom subject to 
controlling forces 


to hydrodynamic or 


and parallel with 


arbitrary origin as 


shown in Fig. 1, the linear and angular momentum equations 


may be written as 


— —_» wed > oe 
F-Slmo)=mi =) LOM Ul 


where 


m = mass of ship 


eit 
rT 


~ = mg 
=ul + vf +t wrk 
u = rate of surging 
v = rate of swaying 


r = rate of heaving 


_» CL), 


linear velocity vector 


= angular velocity vector 


=pi tag¢+rk 
p = rate of roll = ¢ 


q = rate of yaw = 6 


r = rate of pitch = P 


Wes 


(2) 








Kt 
I 


force vector acting on the 
= Xi + Yj+ Zk 


= hydrodynamic force on body 


KK P< 
II 


hydrodynamic force on body 


= hydrodynamic force on body 


>} x 


isco point 


=Ki+M/+ NR 


K = rolling moment about x-axis 


body 

(3) 
along x-axis 
along y-axis 


along z2-axis 


= moment vector acting on body about a 


(4) 


M = pitching moment about the y-axis 


N = yawing moment about the 2-axis 


H = angular moment vector 


-—> — _—> 
=Iypit ag + rk 


a 

= 

aos 
lI 


<wyeanG Zerespectively. 


FIGURE 1. Body axes 
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(5) 


MasSsemement ol inertia about Che 2xis 





incerring Houations (2), (3), (4) and (5) in Eq. (1) 
gives: 
® a2 ° a 
x =on{ ut Eur-yy—%e(B+1 ) + ¥, PE-1) + (Pre g)} 


| 


Y=m\yr yu Pur-JelY + P)+ ka (Q1-P +X, (PPT) (6) 
z= mw purr fu —24(P + $)+ %G(rp-F + (12+ P)} 
fe IP t (L.-Ty) gr +m] U(ur+ Pu- Zu) - %e (vt (u.~Pos)} | 


w= Tee (Ly- le) 1p + Unf te (et Por 1V')- Ag (ur tpv- Fh >) 


= T,7 tC ty- TL P§ +m} %( WtTUr-PW)- ge(a+ Pur- rv} P| 


In Eq. (6) and (7), the left hand side represents the 
forces and moments acting on the ship and the right hand 
Side represents the dynamic responses of the ship, 
respectively. The forces and moments acting on the ship 
depend on the geometry of the hull, the motion and the fluid. 


fmeat 1s, 


7 = : 
F (or /7 ) = {f (properties of ship, properties of 


fluid, properties of ship motion) 


Since the properties of ship and fluid are constant for 
a given ship in a given fluid without excitation forces, the 
forces or moments may be expressed as a function of the 


metion as, 


Fn Ir) = f (Moy ¥ (°, bs, P, 8, ¥, 4, UWP, 6,0 
2 ey P, 2. ‘- é, oe, ele ) 


(8) 


iS 








where xX eho! za represent a system of reference axes 


Oo © 


fixed in the earth and 


oe te oe D G Va = orientation parameters 


U, UL W, P.b a 
CZ, vu, W, ie, Br f 


= motion parameters 


J, J, Ba = control surface parameters 


Considering motion in unrestricted calm water, only three 
degrees of freedom are of concern: yaw, Sway and surge. 
Therefore, the forces and moments are independent of the 
change in orientation with respect to the earth axes. 

Furthermore since the forces and moments produced on the 


Ship as a result of 35 and 35 are normally negligible. 


me eel 


Naomi. my, a, vy, o) (9) 


Since the motion is restrained to the horizontal plane 
oomeene Ssurtace ship in calm water as prescribed in the 
previous section, roll, pitching and heaving are taken as 
Zero, 1.e., p=q=w=0 in Eqs. (6) and (7). Furthermore, most 


ships have Ue = Oecoutnat bas.) (6G) and (7) may be written as, 


iene) ) 


rs 
it 


< 
ii 


Wire ru + XY) (10) 


N = Ig ht m{%(t+ur) | 


ine 





From Eqs. (9) and (10) then, these become: 


X: m(u-Ty er) =f (47, UU, 8) 
vy: m( Utter) = Jy (U,V; ¥, u,r,r, 6) (11) 


N: ee ae a ae = Ey yy L, UV, r, é) 
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Pee eaAMiC STABILITY OF THE SHIP 


Dynamite Stability 1S directly related to the magnitude 
of yaw and sway deviations caused by small initial dis- 
Wawances. In many problems related to dynamic stability 
with small disturbances, the linear analysis of motion is 
quite useful and will be applied in the following to the 


present problem. 


A. LINEAR EQUATION OF SHIP MOTION 

tomsiinplrity the equation of ship motion, the following 
assumptions are made: 

(1) The ship has adequate stability. 

epee ine Water is so calm that the rolling, pitching 
and heaving motions are negligible. Therefore, the motion 
G@umehip May be regared as taking place in the horizontal 
plane only. 

(3) Yaw and sway do not affect the forward speed 
appreciably, and small changes in forward speed do not affect 
yaw and sway motions so that the surge equation is decoupled. 

(4) The ship has one rudder at the stern along the 
centerline and the effect of the propeller is disregarded. 

(5) The forces and moments on hull and rudder are 
expressed by means of velocity and acceleration derivatives. 

In order to linearize the equations of ship motion about 
an initial equilibrium: condition of motion with the previous 
assumptions, only the linear terms in the change of the value 
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mmc Variables from equilibrium, i.e. Au, Av, Ar, Au, Av, 
4r and 4Sare considered. 


Prom £0. (9) then, 


Ke X4Xy au +X, 4V+ XpaVt Xp +¢Xp av t Xpa$ 


epee in initial equilibrium U.= eerie = Us = - =: oe =O , 
AYy=VU-Vo =V 
Ar=Y-% =7 
Peer eC Ul 
AV = oNE on 
sp =¢-h =F 
4$=$-§.=58 


Therefore, 


MN =X 4X OUtX VEX, VAX, U +X, VEX, P+X d (12) 


and likewise, 


Y= Nr Yat yyvr yee tH UtAT tS os, 
N= My tNutl tNy VAN, TtNg ht Ng VN TN GS 


where Xo» i and No are the values of the initial equilibrium 


conditions. 


In the initial equiliburium conditions of straight ahead 
motion at constant speed, there are no forces acting on the 


Ship. Hence 


X= % HN = 9 


dig) 





Furthermore, since the ship is symmetric about the x-z plane 


y= Ya = Xr = Xp= Xp=Nu = Na = 0 
and 


Xs =O 


Therefore, 


Y= X4u + on 
YEMWVE Yt MU tHrr ys (14) 
N =NyU+ Nt Ny Utne + No 8 
Linearizing Eq. (10) gives 
X= ai} ( st AU)= (haar Ontay)-76(n+ar)y 
= m(atl- fav-a44I ) 


This 1S simplified further by neglecting second order terms 
so that 


oe Val = 7) U 


and likewise, 


a m(v+tT Uo + Xe) 


(15) 
N= lr t+ MX (V +1) 
From Eqs. (14) and (15) 
(X,-I)u +Xyru =O 
(16) 


(Yq — r+ Vp VY HV -+ Ome) + Ye d = O 
(Ny- m%e) 4 NV+ (NG Tat + (Ny MG Ue) r + Nh =.0 
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For convenience, Eq. (16) is non-dimensionalized. Every 
term of the first and Second equations for the X-force and 
the Y-force in Eq. (16) has the dimensions of force and the 
third equation for moments has the dimension of a moment. 

To non—-dimensionalize these, the force equations are 

a 
divided through by CL iHu’y and the moment equation is divided 

Sena | os . 
by (2 LHU) where f= dinsity of ship, L = ship's length, 
U = ship speed and H = draft of ship. This gives, 

/ ' an eer if 
if 4 4? fy 4 ey eres i ,o+ / / 
(Yi-m)o +, V+ Mpa )r +(y,-mur+y% § oe) 
1, oe Aliph: ’ eet / a ae , / 
(Ny ~m7G)0 +Ny VY +N; Tp) PON, 7 %G Ue Neg S 
where 
‘ 
y) 4 2 i - 7% UT 
m’zm [LECH ; y= uy io ie ae 


Us = vA, | for small perturbations, 
=] 21H a (fu Ye Ay 

X= n/t Y= MA cnv? , Alf fran? 
Norn M/Lemu, Yea UPEEAY, hm MefZUHY 
Nj=Nr/ E024 U Y= V/2vn Ny mee 
Nj = Ne [Lote 


Zk 





In Eq. (17) the dimensionless time is defined as 


/ EE 
t=7 7% 


B. DYNAMIC STABILITY WITH FIXED CONTROL 
Taking the Laplace Transform of Eq. (17) with the rudder 


fixed at Ne O gives: 


[cY-m)s+ Vi Jutsys [Cy-m %) s +(¥,-m Un’) | r¢s)=0 
ye : (18) 
[Nyon %4)s tive) +[(Nr- J; Js + (Myon 7% u,')] r(s)20 


Then, using Cramer's rule 


O 
) = eS Tena 
v(s A 
nes y= an 
LA\ 
where 
eh : y) / pe / " 
(Ye -) s + My (ye -m % stl Yy-m) 
Xe ee 


(ni - Mi yay aa | (Ny -J,)s (Ny Ge) 


Expanding this determinant gives the characteristic 


equation 


AzAs+BstC (19) 
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where 
a =(e-me- Te )-l Yet Ny %G ) 
B=( Vo - UN; a %) — Yj, ( Ny ~I3 )- 
Ny (V9 Ke) — (Ve a Vig “He? (20) 


nf t 7 / / 
oe VO aa HG )-Ny, (Ye 1) 
The roots of Eq. (19) are 


oo) ae B a 
eee /B-4AC ae (B)-4(&) 
| = a a ad 
2A Z iis 
1B //B tC 
-A VA) 4G) 
2 


Wy == 


Therefore, the free responses of ship, i.e., sway and 


yaw motions are given by 


/ 


2 BUTS Borges Gt 
a € ig 1S (22) 


ae qt at 
aoweane a De 


where v r, and r, are arbitrary constants depending on 


ea 2 
mercial conditions. 

As shown in Eq. (22), the dynamical stability depends on 
the roots of the characteristic equation. Since d7 > 02 ; 


Ain vadj=0 and Lim rit) =o 


L'a t-700 
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aceon. at Gi; <ae 

If (0 then, Vw) + O and r (a) * O which means that the ship 
is dynamically unstable in straight line motion. The root, 
6; is called the ‘stability index' 

In Appendix B and computer program 1, a five foot Series 
60 model of the Davidson Laboratory with block coefficient, 
Cp = 0.7, and no propeller, was modified to be (a) Sep Ie 
(py) Marginally unstable (c) unstable (d) very unstable and 
(e) extremely unstable. These models are referred to as 
me@el A, B, C, D and HE, respectively. 

The hydrodynamic derivatives and stability indices of 
these models are in Table 1. 

Eqs. (21) and (22) were simulated to obtain the response 
of each model to small a disturbance. (See computer program 
foe Fie. 2) 

As shown in Fig. 2 (which is the result of computer 
program 7), the yaw rate, r, as well as the sway velocity, 
v, of the stable model A decreases as time increases but 
the unstable models B, C, D and E tend to diverge due to 


the positive stability index. 


TABLE I. HYDRODYNAMIC DERIVATIVES AND 
STABILITY INDICES OF THE FIVE MODELS. 


acaba cara Ca 


~,0635| —. 280 |+. 0040 |+ eae 2 -.3213 













ef 

038 

Ce fara nas ores + 0b2{ = eee 
es. odib [0h 35| 1268 
oL. | 








71038 
—, (367 |._ A771 > 0065 }— 
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C. STABILIZATION OF THE DYNAMICALLY UNSTABLE SHIP WITH 
FIXED RUDDER 
As described in the previous section, a dynamically 


unstable ship has 


ee a 


B es @ 
rae as =) -4 
(eal ax ero 


To stabilize this unstable ship, @; must be made less 


than or equal to zero, i.e., (<0 or, 


_ B 
KBy.4) oz 


By the nature of hydrodynamic derivatives for the ship 
A>O and B> 0 always. Therefore, C must be greater than or 
emeak to Zero for stability. Thus, C is called the "stability 
ereverion." 


mreom Ha. (20) 


C= (N° %4)-Nol(~ 1) ZO (200) 


Then, considering the characteristics of each of the 
myaroadynamic derivative, mass and the distance of the C.G. 


from the origin of the fixed coordinates on the ship. 


is always large negative quantity 


Me 
N is always large negative quantity 
N may be positive or negative 

Y 


is small but may be positive or negative 
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m is always large and positive 
Xa is always small but may be positive or 
negative 

Paina lly ,; No igs. tne sreatesy Ciliect on dynamical stability 
Since the term (Y-m’) is a large negative quantity in 
eomeation (20a). If, for example, Ny becomes more positive, 
mie stabllity criterion C will become more positive so that 
the stability of the ship will be improved. It is noted 
that Ne can be controlled by a lifting surface. 

Me eociStiuimity Of Dymamic Stability to a Lifting 

surface 
Equation (21) was simulated using the computer 

program DSL/360 Language for the effect of the location, size 
pm@oeaspect ratio of a lifting surface (movable fin) on the 
models A, B, C, D and E. (See the computer program 2). The 
results are summarized in Table II and IIl. 

The following are indicated by Tables II and III and by 
Imeeee 3, 4 and 5: 

ie Lhe position of the fin at which the stability index 
is maximum is independent of the area and aspect ratio of the 
aT), “SORMIRe ETELS POSPLIOn 1S Chie pivot point at which the 
relative flow is parallel to the centerline of the ship. 
(See Fig. 14) 

The parallel flow makes no attack angle on the fin so 
that Yo Nu. Y.. and Ae of the fin are zero but Gee Nae ae and 
i are not zero as described in Appendix A. Therefore, the 


Sooblity criterion C in Eq. (20) is the same as that of the 


2 





TABLE II. EFFECT OF FIN AREA AND FIN LOCATION WITH 
FIXED ASPECT RATIO ON DYNAMICAL STABILITY 


Er 
— 32079] + 0242) ; 
— 0341] + .2026, ; 
(0.34) | (0.305) 

88 

p107 6203 

| (0. 308) 




















pivot 


point 


4- 





foal 











ap al 










7921} -— .6486 | - 
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Meemes the effect of fin location and fin area with 
constant aspect ratio on Stability of the ship. 
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FIGURE 4. The effect of fin location and aspect ratio with 
Conctane area Of fin on the stability of the ship. 
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fo—eiieconditwon. In other words, the stability eriterion, 
Seroeminimum at the pivot point during shifting of the fin 
from stern to bow. 

mits pivot point 1S the position of the fin at which it 
exerts maximum destabilizing influences on Ne ship. 

2) The effect of the fin on stability increases with 
fn Size. 

3) The more unstable the ship, the closer the pivot 
peimt moves toward the midship. 

4) Compared to the no-fin condition, the placement of 
the fin near the pivot point has a negative effect on 
muaprlity. This criterion is called the "destabilizing area” 
in this thesis. This destabilizing area is due to the 
feamenlve change in A, B, and C in Ea. (21) as the fin location 
imei red. C is almost constant near the pivot point and 
the destabilizing area is independent of the size and type 
of the fin. 

5) Moving the fin toward the stern from midship 
improves the dynamic Stability of the unstable ship. Con- 
Wemeely, moving the fin forward of midship does not cause 
much change pont no-fin condition regardless of the size or 
type of the fin. 

6) Stability is sensitive to the fin area but rather 
insensitive to changes in the aspect ratio of the fin. 

Pee ceneotEIVvinty Of Stability of the Unstable Ship to the 

Posmevon. 01 the C,G. ot the Ship 
itetswelear that moving the C.G. toward the bow tends 


to make 'C' more positive in Eq. (20). The effect of moving 


3ke, 





eae was Suudied by use of Eqs. (20) and (21) using DSL/360 
Language for model C. The resultS are shown in Table IV. 


(See computer program 3). 


Pe pie ie tener Or LOCATION OF -C.G. ON 
THE STABILITY INDEX 


Paolo 4022514 .0775 | + .015 | = .0473 |= 2229 


As shown in the above table, stability was improved by 





moving the C.G. toward the bow as expected. However, this 


is not a practical method of stabilizing an unstable ship. 


D. STABILIZATION OF THE DYNAMICALLY UNSTABLE SHIP WITH 

RUDDER CONTROL 

In previous sections, the stabilization of a dynamically 
unstable ship using lifting surface was discussed. However, 
fees clear that the dynamically unstable ship can be made 
directionally stable if the rudder can be forced to compen- 
sate for disturbances so that the ship will maintain course 
with only small oscillations in yaw and sway. If the 
deflection of the rudder, & , is proportional to the heading 
error, Vf, and the angular velocity, r =, then the equation 


of rudder deflection becomes 


Sitv= hbltto+e hPa (28) 
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where k and, CeewcOnsvanes Of proportionality of the 
control system. 
insect tine «iG. 6( 25) amco ba. (17) and taking Laplace 


Transform gives 


lym + vis)+|(Y; ~ nq )s-+ (Yer +% hy) bk] poo 


(24 ) 


lvy-m%.) S +N,| u (5) + [i ~J,) 5 + (N,-m 44th, Ne )s4h No Jabs) = 0 


As these algebraic equations of (24) are homogeneous, a 
mem-trivial solution (Vv 2 + 0) demands the determinant of 
mpiemsystem, i.ée€., characteristic equation, to be zero. 


Let the characteristic equation be 


As+Bs+Cst+D=0 (25) 


where 


A= (Sy-m' Nive Te) MG Ml 7%) 
B=(Ve-m' Ny M2 1K, Me) 4 Yon -E) 
Ny (Yp'- mn %¢) - (Yet h, Vo MG V,~-u%G ) 


Ce (Ye -my(h slg) + Yo Nom +R NE )- 
Ni - 00+ BYE RM NG wi 9G) 


p= bY NS - RM Me 
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/ 7 
Let the roots of Eq. (25) be GG and g;’. Then, 


wit’) 


v(t’) 


/ 


Gt’ Gt’ Gt 
ye +¥%e + Ke’ 


? 


I 


re. Cae Brees (26 ) 
Gt G;, t G3 t 
Ve +U@ +052 


7 


i 


mand only if G0: and ( are nee Av ive> then Yo) and 
Vo) will be zero and the unstable ship will be stable. 


Furthermore, if the roots are three repeated roots, i.e 


Lae } 
é 


Oe Ge = Os, ac; the ship will be stable with the shortest 
possible response. 

By designing for the optimum value of k and in the 
control mechanism, the ship will be automatically controlled 
Porpesuitable stability. 

To find the three repeated negative roots of the 
Characteristic Eq. (25), the 'MITROVIC METHOD' was applied. 
(See Appendix C). 

In computer program 4, the optimum 2, and Rk, for the three 
repeated negative roots of Eq. (25) were calculated for model 
C Bia! FG. (26) was Simulated by the DSL/360 Language. (See 
computer program 4). 

As shown in Fig. 5-1 which is the result of computer 
program 4, the response with proper Rand PF; for three 
negative ont eee Of wWihaEracvTeristle Couation (25) 1s an 
exponential transient. Conversely, if the auto-steering 
device has some #, and /, which give complex roots of charac- 
weristic equation, the response will be stable with some 


Sserliations. 
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The auto-steering device should be able to select the 
Peeper Control costants, P and k, . 
It is evident that the dynamically unstable ship can 


be made to be stable directionally with auto-steering control. 
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IV. MANEUVERABILITY OF THE UNSTABLE SHIP 


To predict the performance of a dynamically unstable ship 
as well as to get realistic predictions of maneuvers with 
large rudder angles, a non linear mathematical model is 


needed. 


A. NON-LINEAR EQUATION OF SHIP MOTION. 

In the non-linear analysis of the Ship motion, the 
following assumptions which are rather typical in the 
literatures Refs. 2 and 8 were applied. 

1) Terms of higher order than third order are neglected. 

2) Rudder force and moment derivatives of higher than 
first order, and effects of rudder angular rate are negligible. 

=) Odd terms in Vv, oe, vy and r in the X-force equation 

and even terms in the same parameters of the Y-force and 
N-moment equations are eliminated due to symmetry. 

4) Cross coupling terms between acceleration and 

velocity are neglected. Taking the Taylor Series expansion 


omeba. (ll) with the above assumptions: 


¥ f(ul. VT, 5) 
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Imeteerating Eq. (27) 
- , 
ult')=f ae’) dt 


/ a / 7 
Cae { ut) dd 
‘ C29)) 


Y(t) =f" Td) dt’ 
with =f rctydt’ 


Note: The time domain should be non-dimensional. 

So far the reference axes x, y and Z are fixed to the 
moving Ship. But it is convenient to take the reference axes 
xX and We fixed to the earth in order to determine the path 


O 
of the ship. Adopting the figure from Ref. 1] gives 
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FIGURE G. Relations between axes fixed 


irae eCarth and axes fixed in the ship. 
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The following relationships may be written as 


X(t) = u(t) Cor WM - vib an vb) 


(30) 
db) = upain YD) +t) cet Pb) 
i 
X.(£) | %,(t) dt 
C31) 
C : - 
Yt) = | 4, h) a 
Non-dimensionalizing Eqs. (30) and (31) using the 
ioral tions, 
Ze Ko . 7 Ky y, Os ee, a LJ 
Ene ; au d 7.q “TW, A =f fe 
eal, | | 
Xo ms Ss Ke U dt 
e.) [_ at = { ip ie 
or 
X(t = [ade . 


Likewise, 


Wet) ={4" dt" 
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pee fURNING CIRCLE WITH SELECTED RUDDER DEFLECTION 
Fig. 7 shows the typical rudder deflection with respect 


to time in a turning maneuver. 
Say) 


oe 


nm Be = ple 


é Ke wom Leng. 
ELGURE “7. hucddes der lecLltion vs. 


time in turning maneuver 








—— — —— — eee =— 


Bsing the ramp function for rudder deflection, then 


4 “1 
S(t) = & i ramp (i Lesa) — (anid cape Mee) (33) 


where 


é = rudder deflection 


Ano Gay) ramp function 
OX = rudder deflection rate 
U tag = time delay 
Xs max time at Evie 
Calculation of steady turning radius. 


The relation of the angular velocity with linear 


velocity gives: 


(radius) x (angular velocity) = (linear velocity) 
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Uf) [war d)+u et) 


or Radius =--— = 


rp) rt) 


Taking non-dimensional terms 


Radue uy furry? _ fy +(wuy 





i U/ 
Chale. Me y() 
Bt / 1 (U')+ (u’)” 
Re. = (34) 
hes 
where / Kadiud. 
Kad = 


Figure 8 shows the typical turning trajectory and 
terminologies corresponding to Table 5. Eqs. (31) and (34) 
were Simulated for the 5 models using computer program 5 in 
turning circle maneuvers with a rudder deflection rate of 
2 degrees per second for the prototype (L = 500') with zero 
time lag and constant final angles of oO 20° and 35° 


Bitemresults are presented in Table 5 and Figs. Sa, 8b and Sc. 
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FIGURE 8a. 





Turning circle of five models with 


35 degrees of rudder angle. 
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curve # : model A 





curve #1: model A, #2 : model B, #3: model 5, #l:rodel D 
ard #5: rodel F 


a) 
(a) 
wD 





x, / 1 


FIGURE 8b. Turning circle of five models with 
25 degrees of rudder angle. 
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curve fl:rmodel A, #2:roJel 8, 43:rodel 3, #u:rojel D, #S5:roiel FE 





FIGURE 8c. Turning circle of five models with 


10 degrees of rudder angle. 
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TABLE V. CHARACTERISTICS OF MODEL A,B,C,D AND E 
IN TURNING CIRCLE 
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340490 
6.9411 
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Note: Adv. = Advance/L 
Tran. = Transfer/L 
T.D. = Tactical Diameter/L 


S.T.R. = Steady turn radius/L 
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Prscussion Of turning circle maneuver 

eens Shown in Fig. 9, the dynamically less stable ship 
has a smaller tactical diameter, advance, steady turning 
radius and circling time than the more stable ship, because 
the yawing rate of the unstable ship is greater than that of 
stable ship at the same rudder deflection. However, the 
turning circle parameters are becoming constant as the 
stability index increases. 

Pee eos Shown in Figs. 10, 11, 12 and 13, the rate of 
change of the turning circle parameters versus the rudder 
deflection is high at the small rudder angle. In other 
words, the above parameters are less sensitive to rudder 
angle as the rudder deflection increases. 

3) The more unstable the ship, the less sensitive the 
above parameters are to the rudder deflection. 

1. The Effect of the Location of the Rudder on the 

Turning Circle 

A change in the location of the rudder has multiple 
effects on the ship motion. One is it's effect on dynamical 
stability, which is the same as the effect of a movable fin 
described in section (II-C-l1) and the other is it's effect 
On maneuverability. 

As described in Appendix A and Ref. 1, 

1) ee nes apo) ) (sake sandependenGlol the=position of 
ime rudder. 

2) N., Yo Ne and -¥ ~ are less positive if the rudder 


meves forward from the stern. 


o4 








2S ) ee and Nia are less negative if the rudder moves from 
the stern to the midship position and less positive from 
midship to the bow. 

4) Since the relative flow passing the ship depends on 
the position considered, the attack angle depends on the 
[Mecarion Of the rudder. Forward of the pivot point the 
attack angle is larger than rudder angle, at the pivot 
point the attack angle is equal to rudder angle and aft of 
the pivot point the attack angle is smaller than the rudder 
angle. 

5) N strongly depends on the location of the rudder 
cance N< = Y¢ -X¢ where Xf Me ene dievance ol the rudder 
position from midship and is negative aft of midship. Ne 
is less negative if the rudder moves forward to the bow, 
but Nod is less positive if 6 is negative. This less 
positive effect decreases the yawing rate, r, due to 
decreases of f., in Eqs. (23) and (28). Decreasing r in 
turn causes the turning radius to increases. 

On the other hand, the movable fin effect of shifting 
the rudder forward makes the ship more unstable which makes 
the ine adius decrease. These two compromising effects 
are shown in Fig. 14 and 15 which were obtained from 
momputer program 6. In computer program 6, model C was 
simulated for the evaluation of the effect of the rudder 
position on turning maneuvers. 


As shown on Fig. 14 which represents the effect of the 


rudder position on the turning radius without replacement of 


OO 











TABLE VI. THE EFFECT OF RUDDER POSITION ON THE TURNING 
CIRCLE OF MODEL C WITHOUT A FIN AT THE STERN. 


os | " 


8 
=—O2165¢1 0.45764 Lege SS) 
-0O.1701 0.4688 2.164 0.59482 
Sie? ao 0.4734 2.144 0.6492 
.64 





-0.1743 0.4739 2.14 
6.64 -0.17556 0.4742 2.1412 O07 056 
6 


-O0.10 aS /_ 
0 6.76 19 6 Soe 
a2 | 







-0.1784 0.424 ee 

2. 
1 | 10 2.268 | 0.00969 
| 

“0.1885 2380 | 0.9301 
-0.18869 2.4 0.9410 
jaaia | 2s | 0.9324 
-0.01396 tech ers’ 0.9445 


a 










-Q.5260 1a 








T/sutper Suyuany Apeeys 














rudder position / L 


the position of the rudier or 


Ta) 
2b 


to 
turning radius without a fin at the stern. 


iy 
w 


Fig.lh Tne effe 


o7¢ 





Tapia la. Tae EPPECT OF RUDDER POSITION ON 
MiG TURNING CIRCLE WITH A STABILIZING FIN AT THE STERN. 














| 
: peed 
Otro LZ ¢ 

023056 ; 

EOL OOr 
0.2810 3.000 
: ; 0.2220 +0.2948 

+0.0874 +O <2520 


Caves 














o8 








Turhines raecaus /2r 














suddenly reversing 
WOun tran. Ul mas 7 
Gabe Oct Fe i | | 
-~-.3 -.2 -.1 O aul ey: .4 28, 
Rudder position / L 


FEGURE lo. The effect of the rudder position on turning 


eirele of model C with a stabilizing fin at the stern. 
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mem tn at the Original rudder location, the steady turning 
radius decreased as the rudder was moved from -0.5L to 
-0.2L, and increased as the rudder was moved from -0.2L to 
meee. AL the rudder position 0.24L, the turning radius 
decreased suddenly from 2.415L to 1.963L. Also, since the 
turning moment at this point was reversed in direction, the 
eaepecsuddenly swings to the opposite side. This phenomenon 
meme Characteristic of the unstable ship as described in 
section IV-D. 

If the turning moment is decreased due to shifting the 
position of the rudder, the yawing rate decreases until the 
ipucider Change occurs. After changing the turning direction, 
the magnitude of yaw rate was suddenly increased. Finally, 
the turning radius was decreased markedly. 

ieestind the Optimum position of the rudder for proper 
stability and maneuverability of model C, the rudder was 
moved forward to the bow and in it's location a stabilizing 


immewas installed. 





As shcwn in we. AEs, and Table Via, the turning radius 
at rudder position -0.33L is a minimum before a sudden 
change of Peete MidwinnaAGdalt POletie ws tani ty Andex 1S 
Meeeacive. Although this minimum turning radius (3.223L) 1s 
ereaver than that of the original model C (2.296L), it is 
meecceptable optimum rudder position due to it’s better 


euabi lity index, i.e., negative index. 
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C. ZIG-ZAG MANEUVER 
According to the standard zig-zag maneuver, the relation 


omene rudder deflection Cy and ship's heading CY) nS) wens 


@escrsbed in Fig. lo. 


Jt) 
J NOAaX 


a A MOK 





FIGURE 16. Relations between rudder deflection, 


ship's heading and time in Standard zig-zag maneuver. 


Using a ramp function for the rudder displacement, the 


rudder angle iS given by: 
SID) =& [romp (tt) - Yau tty) - van? (F-3) + 
Youn? (i-4,) t romp (Fe) | 


- tf se i 
fe i [ Cycke uniet A Ie ) 
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ie = time delay of rudder 


rudder rate (degrec/sec. ) 


ramp 1Unetlon 
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Rudder Deflection and Ship's Heading in 4Z-maneuver. 
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Yawing Rate and Sway Velocity in Z-maneuver 
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FIGURE 19. Path Trajectory in Z-maneuver. 
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In computer program 8, the standard zig-zag maneuvers 
were simulated tor model A and C. As shown in Figs. 18, 19 
and 20 which represent the results of the program 8, the 
following can be observed: 

1) The amplitude of overshoot in yaw angle of the 
unstable ship is larger than that of the stable ship. 

2) Although the first reversing execution of the rudder 
eeamuie Unstable ship is earlier than for the stable ship, 
the period of the stable ship is shorter than that of the 
unstable ship. 

3) Amplitude of overshoot of path in the unstable ship 
fweweoreater than stable ship. 

fee tie reason for 1), 2) and 3) is that the yaw rate of 
the unstalbe ship is greater than that of the stable ship but 
the sway velocities are almost the same for either ship so 
that it takes longer for the unstable ship to get the reverse 


heading after reverSing the rudder angle. 


D. SPIRAL MANEUVER 

Mae standard spiral maneuver of model C was carried out 
using DSL/360 imecomputcer proeran 9. Thre process was 
ferred with 10 degrees port rudder and 1 degree of deflect- 
ing increment and ended with 10 degrees starboard rudder. 
The process was repeated. 

Heoure 20 and Table 7 show the results of computer 
Prorcram 9. 

During the spiral test of model C, the angular velocity 


was reduced suddenly from r' = 0.19 to r' = 0.28 at Be 1a. 


65 


4aw Waves | 





Rudder Angle (Radian) 


FIGURE 20. 


Spiral Maneuver of Model C 
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TABEE VID SPIRAL MANEUVER OF MODEL C 


RUDDER ANGLE YAW RATE 
(DEGREE ) (RAD. /SEC*L/U) 
kok KK KKK 0.36548072 
- 9.000 0.35614842 
- 8.000 0.34629613 
-~ 7.000 OPSsooDlO | 
- 6.000 02524 46226 
- 5.000 0.31269461 
- 4.000 0.29961270 
- 3.000 OF 285135420 
- 2.000 0.26874071 
~ 1.000 O0.24949914 
~ 0.000 O22 5505094 

1.000 0.18941998 
2.000 —-0.28294164 
3.000 ~0.29785007 
4.000 -0.31125194 
5.000 -0.32350868 
6.000 -0.33485866 
7.000 -0.34546858 
8.000 -0.35545886 
9.000 -0.36492395 
10.000 ~0.37393415 
11.000 —~0.38254720 
12.000 -~O0. 39080832 
JONG (S) 20 SOoOOLs L 
10.000 -~0.37394869 
9.000 -0. 36493850 
8.000 -0. 35547382 
7.000 -~0. 34548380 
6.000 —-0.3343 7481 
D> 000 =Qys2 352056 

- 4,000 =—Oe e027 167 
3B. 000 -0. 29787230 
2.000 -0.28296709 
000 -0.26596129 

=" 0..000 -0.24573880 
- 1.000 -0.21966898 
- 2.000 =O) A GV eaves rete: 
- 3.000 OFZ Gol0946 
- 4.000 O29 959 112 
- 5.000 OR ozo 7542 
- 6.000 0.32468486 
- 7.000 Qi She os oeyeee 
~ 8.000 0.34628141 
- 9.000 OeS50 le 546 
So Ba Es 0236026072 
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degrees port side. That means any increase of the rudder angle 
tO port by 1.5 degrees will cause the model C to suddenly 

Swing from starboard to port until r' = 0.28. Likewise, when 
the rudder angle increase from port to starboard the model C 
can turn against its rudder by d= 2.09 degrees Starboard and 
then suddeniy swing in the opposite direction to a new stable 


mesition r’' = 0.28. 
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V. CONCLUSION 


This thesis has presented the characteristics of the 
dynamically unstable ship in stability and maneuverability. 

As discussed in Part III, the most powerful factor in 
Sieapllizing a unstable Ship 1s a lifting surface installed aft 
Sftemidship: the addition of area of fin is more effective 
free tabllity thaw the aSpect ratio. The effect of a lifting 
surface may be obtained by designing a fine stern with neat 
flow lines and deadwood or a proper rudder. It was found 
that the fin installed near the pivot point destabilizes 
the ship. 

As discussed in part IV, the unstable ship is more 
maneuverable in a turning maneuver but less maneuverable in 
the standard zig-zag maneuver than the stable ship. In 
turning maneuver, the more unstable Ship was less sensitive 
to the rudder deflection. 

As discussed in section iil Ne followings are 
considered necessary for optimum stability and maneuverability 
ea Ship. 

IDmeuecalting the rudder at the point aft of the midship 
where the steady turning radius is minimum. 

2) Installing a stabilizing fin or deadwood as the 
stern in order to maintain adequate stability. 

DSL/360 Language was found to be very powerful for 


Simulating ship motion. 


8, 





VI. RECOMMENDATION 


All results in this thesis were obtained by simulation 
moe DSL/360 without real model tests. 

The guthor would like to compare these results from 
real model tests in the future. The results in the thesis 
are in non-dimensional values so that they can be converted 
to full scale ships or models. 

Also the following itmes are recommended for future study; 

1) Maneuverability of the unstable Ship in Z-maneuvers 
other than the standard Z-maneuver. 

2) Maneuverability of a dynamically unstable ship with 
auto-steering control. 

3) Ship design for the optimum stability and maneuver- 


ability. 
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APPENDIX A 


Prediction of Hydrodynamic Coefficients in Linearized Equa- 


jmen of the Ship. 


Using the method in Ref. 1 and 2, the following predic- 


tions can be obtained. 


(1) A fin travelling at a forward velocity u and trans- 


verse velocity v has 


(Yr), = Ag (2), 


ae (34) 
where 
Ag = “4 /LH 
At= Qa of fw 
Goa it Cooffeccont 
A= buff ang. 


oO CL = as 


i 





~~ 


2/3 oe for general aspect ratio fins 
tag 


= > Ay for very low eS jie Ny ya ro at SS 
(Mp = Xe (My He 

OM Ve = Xe (My De (34a) 
(Aye = (4 OW) ¢ 


V1 





(Nee = eae ual 








where 
CL = aspect ratio 


Coy = effective aspect ratio 
DL 
om ost 
hy 


hi = atihecwiwe Span 


(Mis = Xe (Yi Ie (35a) 
” a vad oe 
(YP Ie a ip Vu Me 
Je ae 2 r 
(Ne He = (KE) (% )- 
Considering the vudder as a fin, 


= “OM Is 
Ne =~ 3 Ve? 


(35b) 


V2 





(2) Considering the ship's hull as a fin with consider- 


able drag coefficient C 


(Kid, =—[ roe Or (36) 


where ITH 


w | 
=| 
| 
N| 
= 
nS a 
| 
- 


(Ny), = -Cin- R, 1')* He Vos Mp 
( yr I oe -R m + Xp Yo Dh 
(Nyy = - Mg Ve D5 (Vy ), 


(36a) 
Po Sih, 
d Stent. p : 


bow 





car 


(Y¥;)) eS IO 


(see Ref. 1 for more details) 


(6: 





(3) Assuming a large deadwood which has a sufficiently 
Meveaspeck: ratio, an effective span equal to the draft of the 


Ship and a triangle shape, than 


H 


(Yor Der — — 
CMy ler = - Eye )¢ 
OM er =- EON VE 
CMe dg = Yr) 


CY.) £1 = mAs he 
(Dex “ag I) yA 


(37) 


(N., Me = =-t(y, de 
C fie =~ 2 (Yip) pi 


CME) Er = al Yo Me 
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(4) Assuming that a rudder or small deadwood has 








BCL _ 21 
ae - I+ 4z (ie lOosunc ii measpece ratio, then 
“of 


, Lit 
My ey = Age ( ) + deg ) 


(Nv dyn = "2 (Wen 


Wa =-2 (Yr) 2 
(Ny )f2 = (Ny deo 


(38) 


( Vi) —— AN At2 he2 
v i (Glan sei) 


i 


(My )e2 a (Yo Mp2 
C¥elg2 =-F(%y )p2 
(Np lp = 4 (Yo Mes 


I 


(see Ref. 1 and 2 for more details) 





APPENDIX B 


Models for Stable, Marginally Unstable, Unstable, very 


Unstable and Extremely Unstable Ship. 


In this section it is intended to obtain models for 
several types of unstable ships. These will be simulated 
Weine DSL/260 to study their stability and maneuverability. 
The original model was adopted from a five foot Series 60 
model (with block coefficient 0.7 and no propeller) as 


originally developed by the Davidson Laboratory. 


Principal characteristics of parent model. 


Lee eae |? a Sa 2 Peer, a os Pr en & as a ae 
ESTOS GUE) O sl CNT ke, ae eae se Oia ea 6 Ba a 2 ieee a 
JUNE SESS GIS) ts een 9 ee ee On 26 (—h 1h. 
Block coefficient (C,) Ree Rr  Or rt a net 0.7 

LCG from bow (X,) ia eee assine, Se Wee es. Sep lS veri es 
Displacement (A) REE ae hy NE gia che en ae 41.64 lbs. 
Area of rudder (Ap) oe ap ee ee O202 Sais ik: 
Rudder span Ee ee eta aie 45 ofp a G22, 22 

27UuC (CIENC —C1d C166 Ree aera ee aa OAS 70 es i a7 
WATE SHS) CTT) 5 Aa a 1.292 slugs 
RUBS. (CLO TESIOTE C112 Van nent ep na O...2 


Longitudinal added mass 
CoOcit . (kjm' ) RM Es ooo, a5 Bd dc 20s aw we 0.004 


Pabeta ieee cminasoecOectt. (M' )s....... CxS 
Ph 
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Rotational added mass coeff. oe oe OZNG5S 
ECE) Ot fo eietel yakropay 10) see) aba) eee OezZo] 
es Gy otistsareral added mass from 


TE 
midship Gara A coal” eg A tt a 0.024 


Diet COecrt. at Zero drift angle ....... 0.019 


By using Appendix A and computer program 1, the parent 
model can be modified to the desired models. From the 
results of computer program 1, the five models were selected, 
which were named stable model A, marginally unstable model B, 
unstable model C, very unstable model D and extremely un- 
stable model E. (See computer program 1, Appendix A and 
Table I) 

The hydrodynamic derivatives in Table I were considered 


to be inherent for each model in this thesis. 
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AP VEN DEX Cc 


pevectron Of vthe Optimum Values of the Proportionalities 
of the Rudder Deflection to the Ship's Heading error and 


Mawine Velocity. 


Let the rudder deflection equation be 


SH= Rv + ay, (see Eq. 23) 


The ship motion under such rudder control has the 


characteristic equation, 


3 y 
ee ea ea ta) tee 
JS = ie 4 Cs / Ip. O (see Eq. 25) 








Let 

ee ee 

ae A / ,» ven 

Bo Be eels 
a + —— 2 ne ay pee  ) eee 

ad os ie ae ne 

3 2 a, 
ov Ps Ppt ge =o 


OL (Da po 5, P T b, =) (39) 
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where 


gee Ac AD 
Tee and a= Be (40) 


From Eqs. (24) and (25) 


4 ie ene f / 4 / ? , 
A = (jm Nhe BI-Up me Mihm) 


(41) 
Z = G ke, + Cs 
where | 
C= Mg (I )-Ye (We '%_) 
/ y ; y, / f ; ‘ a / 
¢, = (%y male -an He) + Voy (oly Bede (yp tn) 
Likewise 
C2 Gk te K+ Cy (42) 
where 


cy = Wig (Yg.-m') Ye Me) 
ae Ny Ve Ng 


/ / / / / / i, 
Cys My (Mp - He )- Mp Yr le) 


Dez Oy R, (43) 
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where 
7, / 
C= (Or Ns ~My Yo) 
iso. (41) gives 
2 2 vi 2. 
BB atom ee, Ge ee és 


2 -: 


a 
ae 00. See = 4 . 
ie © R> -k oC, Ck, $360, RT Ce 


Eq. (40) gives 


aa) 
2. ” yy } ae ‘ 
b, Ge “2, 2G; 2 Ao Ks } saz Ee Re oi C4 K+ Cs ) 


Eq. (40) gives 


> eee 
pe ib 
ee 
was ae p> 
— 


2 


3p 3 MI ae ee | > ae, 
bp (CG R, pee CG cont DC) Co h, 0, )=A CER, 


Thus 


3 3 2: 2) 2p 3 
: eee ys } i. 
R, = ae et 2 Giga? CC Rio ) 


my (45) 
Igpeertine Bq. (45) into Eq. (44) 
a 3 > 
BR, +B; b+ Bi h+ B, = O 
where 
3 
By, = G C3 Do / ef 
Bie) Gye. D) (46) 
Be = cae os b, C* 


CG, A 
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50, be Do. 1 Glee 
n= ea lf i 





ee 
i, =: G2. G Oo +ACr-b, C 


fMaebq: (46) all coefficients except bo eval by can be 
obtained from hydrod-namic derivatives of a given ship. 


Saiculator of ie and by 


For 3 real negative roots of the characterisitic 


equator, As? at Bs“ tees tee tne point M on Mitrovics 


chart (see Ref. 7) should be in $= 1.0 region 





for 3 repeated roots, M should be the cusp point. 


ME tne Cusp point , 


dl be 
a5 


Mbp Loof dine 


a 


dy dh | du. 


= undefined 


Therefore 


df 0 db, _ 
Ae =a) and re 


ok 





ToOrutteombiicamerder characteristic Eq. with = 1.0 


be = PD. we 
b, = Z2U)4 -3 we 


ees = 0 
of Ne. 


[ 
O64 Utz= Zz, 
| =e 
O05 b, = eh and b, a 2 


(see Ref. 7 for more details) 


Ure t ew LuslSedliirenult to find the exact cusp point by 
the digital computer due to the non-integer 


values of De and b,- 
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< COMPUTER PROGRAM #2 > 
mk EVALUATION GF THE EFFECT OF FIN LOCATION AND SIZE ##% 
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